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During last fifteen years one can observe explosively 
growing interest in rigorous theoretical results of non- 
equilibrium statistical physics implied by time reversibil- 
ity of microscopic motion. In part this is caused by ap- 
pearance of new possibilities for justification of the theory 
on mesoscopic level. Several such experiments and some 
aspects of the underlying theory were reviewed one year 
ago in Physics-Uspekhi l]. However, a state of the the- 
ory as the whole was reflected by [l| not quite correctly. 
The matter is as follows. 

Interest in the mentioned rigorous results has very 
long history. It is sufficient to recollect the Kirch- 
hoff law [if, Einstein relation 0, 0], Nyquist formula 
[H |(| and the unifying fluctuation-dissipation theorem 
(FDT) @, 0. Later, the Efremov's "quadratic FDT" 
8] and Stratonovich's "four-index relations" [9, 10] had 
appeared. We have displayed our interest in 1977-1981 



in series of works where for the first time obtained [11 1 
and investigated [l2T - fil| the "generalized fluctuation- 
dissipation relations" (FDR), or theorems [Hj], which 
universally connect probabilities of observation of mutu- 
ally time-reversed processes and entropy increments of a 
concerned system in these processes. The first of such re- 
lations was formula (7) from [ll |. In the same works and 
in [l6l - |22j we, and then in [23l428| one of us, considered 
also various modifications and applications of FDR. 

Undoubtedly, many new representations and applica- 
tions of FDR may arise with time. However, hardly one 
can rise the degree of generality of FDR as written in 
their abstract form, - e.g. like equality (2) from [l7| . 



p(n+) ex P [-A5(n+)] = p(n_) 



(i) 
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(see paragraph III.B), - since such relations already ab- 
sorbed all given by time reversibility of Hamiltonian 
micro-dynamics. Therefore, similar theoretical results of 
authors of the today's "boom" , - started in 1997-1999 
from works by Jarzynski [3ll - [33j and Crooks f34U36j |. - 
add nothing to our FDR in the sense of generality. 

Nevertheless, in current literature, even including re- 
views [13] (see also references in [29|, H3, [H, [39|), the 
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opinion propagates as if our results are particular cases of 
achievements of "new age" . At that, as a substantiation 
of this opinion, it is affirmed that we resorted to restric- 
tive conditions, namely, coincidence of time-dependent 
Hamiltonian parameters at beginning and end of a time 
interval under consideration. Though it is easy to certain 
of falsity of this affirmation, merely thumb ing through 
our papers (especially since main of them - |lll4l3| | - are 
freely available from internet). To be precise, in review 
[HI and paper [39| (see also references therein) a differ- 
ent, rather objective, look at history of the subject is 
expressed. But it is not dominating, and, unfortunately, 
the misunderstanding had penetrated into review |l| too. 

At this point, it is necessary to underline that our for- 
mulae and Jarzynski and Crooks formulae have different 
statistical meanings. At that, the first directly embrace 
not only closed systems, - to which in fact the second 
are addressed, - but also open systems. In spite of this 
essential difference, the first and the second outwardly re- 
semble one another. Apparently, just their resemblance 
causes inadequate perception of our formulae by some au- 
thors and thus the mentioned misunderstandings. There- 
fore, their "victims" can see only a part of true contents 
of the theory. As the consequence, they do not see many 
applications of the theory and ways of its justification in 
experiments like that described in 

These circumstances induced us to write the present 
"methodical remarks". Below, we shall try to prove, in 
visual fashion, all the aforesaid about relationship of our 
results and results of Jarzynski and Crooks (including 
the so-called "fluctuation theorems" ) and about origin of 
the "misunderstandings" . Simultaneously, we shall point 
out key aspects of the generalized FDR, not pretending 
to review of all of their forms, but illustrating their pos- 
sibilities by selected examples. 

We hope that the material we suggest will be useful for 
all interested readers. For simplicity and brevity, we con- 
fine our consideration to classical statistical mechanics, 
partly following our critical notes [29|, [3(| and the review 
[l| which stimulated our activity (we are grateful to its 
author therefor). Specific remarks on quantum case can 
be found in [30j . 



II. OBSERVATIONS OF NON-EQUILIBRIUM 
SYSTEMS, THE LIOUVILLE THEOREM, AND 
STATISTICAL EQUALITIES 

A. Objects under consideration 

The matter concerns thermodynamically non- 
equilibrium dynamical systems subject to external 
influences. Speaking about evolution of such system 
from viewpoint of statistical mechanics, anyway one can 
not do without Hamiltonian of a system, H(q,p;t), 
and probability distribution of its microscopic states, 
D(q,p;t) , where {q,p} = T are system's canonical 
variables. Time dependence of Hamiltonian can be 



written also as H(q,p,x(i)) where x(t) represent some 
external forces, or fields, or potentials, or governing 
parameters. At that, nowadays it is popular to term 
time dependence of x(t) "protocol". 

Usually it is assumed that if a system stays left to his 
own for sufficiently long time, - i.e. external forces keep 
constant, x(t) = xq = const , and do not change its en- 
ergy, - then one can prescribe for it the Gibbs thermody- 
namically equilibrium probability distribution function, 

D eq (T,x ) = exp{[F(z ) - H(T,x )]/T} , (2) 

where T is resulting temperature of the system, and 
F(xq) its free energy determined by the normalization 
condition: j D(T;t)dr = I. Obviously, this "axiom" 
presumes that the system either contains a large enough 
(ideally, infinite) "thermostat" or was contacted with 
such one. 



B. The Jarzynski equality (JE) and related 
misunderstanding 

Let us address to "misunderstandings" , pointed out in 
Introduction, as they have printed in [lj . More precisely, 
to comparison in []J between "Jarzynski equality" from 
[3ll | (formula (9) in [y) and "Bochkov-Kuzovlev equality" 
(evidently, formula (8) from [iT|). 

In [3lj . a system was considered which initially, at 
time moment t = 0, is in equilibrium state © and then 
evolves being led out of equilibrium by changing of ex- 
ternal forces from x(0) — xq to x{9) at some moment 
9 > when observation of the system is finished, (in [3l[ 
the forces are designated as A(t) while 9 as tf). The 
Jarzynski equality from [31| reads 

( exp (- W/T) ) = exp (-AF/T) , (3) 
W = H(T(6),x(e))-H(T,x(0)) , (4) 

where W is change of full energy of the system during 
observation time interval; AF = F(x(9))-F(x(0)) ; T(t) 
is current system's microstate treated as strictly deter- 
ministic function of its initial microstate T = T(0) and 
functional of preceding forces' values; the angle brack- 
ets denote averaging over given distribution of initial mi- 
crostates, (...)= f ... D(T; 0) dT , in the present case - 
over D eq (T,x(0)). 

In this respect, in [l| it is said: "If final value of pa- 
rameter A coincides with initial one, A(i/) = A(0), i.e. 
the process is cyclic, then AF = , and equality (9) re- 
duces to that earlier obtained by Bochkov and Kuzovlev: 
(exp(-/3W0) =1". 

From this statement any reader must conclude, of 
course, that the Jarzynski equality is more general and, 
therefore, study of the theme should begin from [3lj . 

However, if author of [l| reproduced derivations of both 
the equalities, not one of them only, then it would be 
clear that, firstly, they concern different observations 
(quantities) and even, generally speaking, different initial 
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distributions, so that their literal comparison (by their 
surface appearance) has no sense. Secondly, in [ll| (as 
well as in other our works) thete is nothing similar to the 
condition A(tj) = A(0) (moreover, time dependence of 
forces there has no restrictions at all). 

Let us consider the actual state of affairs more care- 
fully. 



C. The Bochkov-Kuzovlev equality (BKE) 

We have to emphasize two essential differences of our 
equality ((8) In [Tljj ) from the Jarzynski equality (JE). 

First, in our equality the place of the quantity W is 
occupied by 

E = H (T(d)) - flo(T) , (5) 

while place of initial distribution ([2]) by 

D(T;0) = D m (T) = exp{[F - H (T)}/T} , (6) 

where time-independent function H (T) is formally ar- 
bitrary (within reasonable limits, of course), and Fq is 
determined by normalization (in factor exp [Fq /T] 
was written as l/g). Second, our theory permits not 
only smooth time variations of the forces but also dis- 
continuous variations ("jumps"). 
Our equality reads 

(exp(-E/T)) = 1 (7) 

It can be proved with the help of obvious identity 

exp (-E/T) D m (T) = D m {T{6)) (8) 

and the Liouville theorem, according to which Jacobian 
of micro- variables transformation from T to T(t) al- 
ways equals to unit, dT/dT(t) — 1 (phase volume is 
conserved). Namely, 



-E/Ts 



J e~ E/T D m {T)dT 



A«(r(t)) 



dv 

dT{t) 



dT(t) 
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(since distribution Di n (T) is normalized to unit). 

Clearly, this derivation of BKE is completely indiffer- 
ent to how the forces x(t) behave with time. Replacing 
in these "calculations" — E with AF — W , and D in (T) 
with D eq (T, x(0)) , one gets proof of JE ©. 

In PJJ) the function Hq(T) was chosen in the form 
Hq(T) — H(r,xo) with definite fixed value xq (for in- 
stance, xq = 0), i.e. represents some essential part of 
full Hamiltonian and total energy of the system. In such 
case, the distribution ^ (generally non-equilibrium) co- 
incides with equilibrium distribution @. 

However, in our consideration, in opposite to [31], not 
necessarily x(0) = Xq , that is forces can commit "jumps" 
in the beginning of observation, i.e. can be sharply 
"switched on" (as is natural to say at xq = 0). 



D. Comparison between JE and BKE 

In view of the mentioned possible arbitrariness in 
choice of H (V) , BKE (J7J) can not be deduced from JE 
([3]), while the opposite is easy realizable, just due to the 
arbitrariness [29(. Thus, with regard to formal generality 
the priority belongs to BKE. Nevertheless, after a choice 
of -ffo(r) is made, H (T) — H(T,xq), these equalities 
become independent. It is useful to compare them at 
x(0) = xq , when tboth hey correspond to exactly same 
protocols and statistical ensembles. 

Let us write system's Hamiltonian in the form 



H(T,x)=H (T)-h(r,x) 
h(T,x) = H(T,x a ) -H(T,x) 



(9) 
(10) 



interpreting -ffo(r) as "internal energy" of the system, 
while — h(T, x) as its "energy of interaction" with sources 
of the forces. Then, according to system's equations of 
motion (Hamilton equations), 



W 



dx(t) dh(T(t),x(t)) 



dt 



dx(t) 



E = 



dT(t) dh(T(t),x(t)) 



dt , 



dt , 



dt dT(t) 
W-E = H{T{0),x{6)) - H(T(9), x ) = 
= -h(T(9),x(9)) 



(11) 
(12) 

(13) 



It is clear from these expressions that the difference be- 
tween W and E disappears only in special case x(8) = 
xq , and then JE and BKE "reduce" one to another. But 
in general none of them reduces to other. 

It is seen also that E can be interpreted as energy 
dissipated by the system, and this interpretation of E is 
valid, along formulae (fT2)) and (fT3)) . also at x(0) ^ xq 
(and in presence of any other jumps of forces). 



E. BKE for open systems, or when JE is out of 
work 



What of the two equalities, JE and BKE, is more use- 
ful? The answer depends on character of action of the 
forces onto system. In reality, one of forces do work at 
system only when they are changing with time, while an- 
other can do even when being constant (if only they are 
not "equal to zero"). The JE is more convenient in the 
first case, which demands forces to vary smoothly in the- 
ory, since in practice their sharp variations may be non- 
realizable. Example of such force (parameter) is given 
e.g. by position of the piston at Fig.l in [l[ (though, 
its pulling out can be arbitrary fast). What is for the 
BKE, it is more convenient in the second case, when the- 
oretical jumps ("switchings on" and "switchings off") of 
forces are quite adequate to practical actions. Example 
of such force is given by the current-induced torque at 
Fig.3inQ. 
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When considering the second case, it is natural to es- 
tablish reference point of forces in (fit)]) to be x = , 
thus presuming (without loss of generality) that just at 
x = our system is free from external influences. For 
example, if x(t) represents potential drop through an 
electric conductor, or electric force (field) acting onto a 
particle (charge carrier) in an unbounded (in the "ther- 
modynamic limit") medium. 

At i^O, such a system comes to (quasi-) stationary 
non-equilibrium state with constant dissipation of energy. 
In other words, such system (or some its part) behaves 
as non-equilibrium "open system" . 

Now, we would like to demonstrate that JE becomes 
inadequate in such situations. Let the force be sharply 
"switched on" at initial time moment, i.e. x(t) = for 
t < and x(t) = x = const ^ for t > 0. Then, as 
one can see, the quantity W (jTTJ) as a function of 8 
degenerates into constant W — H(T, x) — H (T 1 0) , since 
at t > total system's energy stays constant. 

Consequently, JE ([3]) gives no information about what 
happens during non-equilibrium processes at t > . 
Moreover, JE loses significance at all, since the average 
there, (exp (-W/T)) = J exp {[F(0) - H(T, x)]/T} dT , 
turns to infinity. For instance, if denoting with Q 
( Q € r ) coordinate of particle in the second of men- 
tioned examples, we have H(T,x) = Hq(T) — xQ . This 
evidently implies exponential divergency of integration 
over Q , so that AF = — oo at x ^ . 

At the same time, E = h(T(6),x)-h(T,x) = x[Q(9)- 
Q(0)] gives energy transferred during all the observa- 
tion time from external force's source to system's inter- 
nal degrees of freedom (to thermostat), i.e. dissipated en- 
ergy. Therefore, the BKE ([7]), along with other FDR (see 
Sec. Ill), brings [T^ - [T5| useful information about interre- 
lations of dissipative and fluctuational characteristics of 
non-equilibrium open system (while its free energy falls 
out of play). 

In this kind of applications it is reasonable, if not nec- 
essary, to specify the interaction Hamiltonian (ITU1) in © 
as bilinear form h(T, x) = x ■ Q(T) , that is write 

H(T,x) = H (T)-x-Q(T) , (14) 

where Q(T) are internal variables "conjugated" with ex- 
ternal forces. Then the dissipated energy (Tl2l) acquires 
visual (and customary) expression 

E= [ I(t) -x[t)dt , (15) 
Jo 

where I(t) — dQ(T(t))/dt are "currents", or "flows", 
conjugated with forces. 



were aware of our works [ll| or [14| then they would 
have motives to verify, in addition to JE ([3]), also BKE 

([7]) , by measuring not only integral W = — J Q M Q dt 

but also integral E = J Q M Qdt (dissipated energy). 
Simultaneously, it would be possible to verify equality 

P(E;x) exp(-E/T) = P(-E;x) , (16) 

in which P(E; x) is probability density distribution of E 
under given protocol x — x(t) , and x(t) — x(6 — t) . This 
is direct consequence of FDR for probability functionals 
(see below p.III.B and p.III.F), with x{t) = M{t) (M is 
the torque) and Q = 9 (angle of rotation). Of course, 
the equality (1 16[) as well is valid without any restrictions 
on time dependence of forces (i.e. torque here). 

By the way, notice that experimental "verification" of 
formally exact inner facts of theory means, in essence, 
that of either paradigm of the theory as the whole (statis- 
tical mechanics) or correctness of posing and description 
("stochastic model') of the experiment itself. 

G. Exact equalities for evolution from 
non-equilibrium states 

Let us consider a system with Hamiltonian of type (TTJJ 
and statistical ensemble determined by non-equilibrium 
(normalized) initial distribution as follows, 

D(T; 0) = exp {/3F(X) - \J3H (T) + X ■ Q(T)}} , (17) 

where (3 — 1/T, Q(T) is some set of individual or col- 
lective variables of the system, and X is conjugated 
set of its internal "thermodynamic forces" character- 
izing its nonequilibrium. Such the distribution corre- 
sponds to maximum informational entropy of the sys- 
tem under given mean values of Q(T) and its inter- 
nal energy -ffo(r). Application of the abstract equality 
(UJ to this situation, - by replacing H (T) in (J7J with 
H (T)+TX-Q(T),- yields 

(exp (-AS)) = 1 , (18) 
AS = (3E + X ■ [Q{T{6)) - Q(T)] = 

= / I{t) ■ \fix{t) + X] dt (19) 
Jo 

In concrete examples the quantity AS usually can be 
treated as increment of entropy of the system during its 
evolution caused by both the initial nonequilibrium and 
external influences. 

In (l3l - tl6| it was shown that FDR corresponding to 
such ensembles create reliable ground for non-linear non- 
equilibrium thermodynamics. 



F. On experimental verification of exact theoretical 
equalities 

Here, we cannot but notice that if authors of the con- 
cerned in [l[ experiments with torsion pendulum (40j 



H. Thermodynamical inequalities 

It is well known that for any random quantity A the 
inequality ( exp A ) > exp (^4) is valid, which turns to 
equality only for A — const . 
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Substituting — E in place of A , with E from ((5j) or 
(|12p or (j!5[) . and combining this inequality with equality 
((TJ) , it is not hard to conclude that (E) > . Hence, if a 
system is initially equilibrium, then on average (ensemble 
average) it always takes from sources of external forces 
and dissipates a positive amount of energy. 

Now, let our system is non-equilibrium already in the 
beginning of observation, like in the preceding paragraph. 
Then equality (fl~8| implies inequality (AS) > . It al- 
lows negative values of the average (E) , that is now the 
system can do (useful) work against external forces, when 
— (E) > Q . At that, the inequality (AS") > together 
with (fT9|) establishes restriction on possible value of this 
work, dependently on the initial system's nonequilibrium. 

Next, we shall go to more informative exact statistical 
equalities including time-reversed process, and consider 



mean value in ensemble of trajectories generated by ini- 
tial distribution ©: 



how results of our works 
particular, by Crooks [3 



HI EH El relate to results (in 
36f) what appeared after 1997. 



III. TIME REVERSIBILITY OF MICROSCOPIC 
MOTION AND GENERALIZED 
FLUCTUATION-DISSIPATION RELATIONS 
(FDR) 

A. Rules of time reversal 

For beginning, let us introduce designation T = 
{?! — p} and assume that H(T, x) — H(T, ex) , where 
for some of forces ("even") e = 1 while for all others 
("odd") e = — 1 . In fact, it always is so, if one suitably 
defines the forces x and includes into their list all static 
odd parameters of Hamiltonian when such parameters 
exist (the standard example is static external magnetic 
field). 

The time reversibility means that any system's phase 
trajectory can be passed back if making at some its "end" 
time point, 9 , replacement of T(9) by T = T(9) and 
x(t) by x(t) = ex(9 — t) . Here T is initial microstate of 
backward motion, and t is new time count off back from 
the old end point. 

At that, in the reversed time, the reversed trajectory 
r(i) = T(9 — t) depends on V and x(t) in exactly the 
same way as original "direct" trajectory T(t) depends on 
r and x(t) . Correspondingly, observation of any vari- 
able Q(T) with definite time parity, - Q(T) = eQ(T) 
(e = ±1), - instead of Q(t) — Q(T(t)) on direct trajec- 
tory will display Q(t) = Q(T(9 -t)) = eQ(9-t) on the 
reversed trajectory. For the "currents" I(t) = dQ(t)/dt 
then one will see I(t) = —el(9 — t) instead of I(t) . 



B. Generalized FDR 

Let ${r(t)} be some (arbitrary) functional of system's 
phase trajectory on time interval (0, 9) . Consider its 



(*{r(t)}) x(t) = J ®{T(t)}D in {T)dT , 



(20) 



where the comment under right angle bracket reminds 
of external conditions influencing the system's evo- 
lution. Along with (|20)) consider ensemble average 
(<!>{r(£)} exp(—E/T)) x u\ . In the phase space integral 
in this average let us (i) use, like above in paragraph 
II. C, the identity ([8]), (ii) go from integration variables 
r to T(9) , applying the Liouville theorem, then (iii) go 
from r(0) to T = T(9) , and finally (iv) apply the above 
"rules" . As the result, we come to equality 

(*{r(t)} exp (-E/T)) x[t) = (*{f (t)}) 5(t) (21) 

Notice that due to the operation (iii) on the right here the 
averaging is performed over initial distribution Di n (T) = 
D eq (T,xo) = D eq (T,xo) (xq = £Xq) which essentially 
(and at xq — or e = 1 literally) does not differs from 
distribution D in (T) on the left. In order to underline 
this principal circumstance, the left and right protocols 
in (f2"Tj) can be written as {xo,x(t)} and {xo,x(t)} (thus 
indicating also preparation of initial states). Below in 
this paragraph we suppose, as in [TTJ] , that xq = . 

Choosing ${r(£)} properly, it is easy to transform 
(TUT) into various "generalized FDR" for characteristic or 
probabilistic functionals. For instance, for some set of 
variables Q(T(t)) let 



*{r(*)}= n s m)-Q(r(t))} = 

P{Q(t);x(t),T} (22) 



o<t<e 



This is probabilistic functional (probability density dis- 
tribution in functional space) of trajectories Q(T(t)) 
which corresponds to concrete initial microstate T . Then 
the average ([20|) , 



(23) 



P[Q;x} = J P{Q;x,T}D in (T) dT 



is probabilistic functional (PF) of trajectories Q(T(t)) 
in the framework of given statistical ensemble of phase 
trajectories. 

Further, let the set Q(T) is wide enough to express in 
its terms the functional E : E = E[Q; x] . Then insertion 
of dm into ([21]), subject to pg|. leads to formula (7) 
from 11], 



P[Q;x] exp(-£J/T) = P[Q;x\ 



(24) 



(in contrast to [11( , here the parities e , e are included by 
the symbol "tilde" ) . Both sides of (|2"4")l relate to one and 
the same (equilibrium) initial distribution D eq (T,0) . 

In case of "bilinear" Hamiltonians (fT4")) the set Q what 
appears in (|24p should cover the set appearing in (|14l) . 
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Therefore the second can be identified with the first if 
setting all additional conjugated forces to be identically 
equal to zero. 

Notice, besides, that at ${r(i)} = exp {-[IF - E]/T} 
and x(0) — xo relation (l2"Tj) turns to the JE ([3j> . 

In order to generalize equalities (|2~T]) and (|24| to 
the "initially non-equilibrium" ensemble introduced in 
p.II.G, it is sufficient to replace E by AS and add 
the thermodynamic forces X to comments under angle 
brackets. The resulting analogue of (IM1) can be repre- 
sented by formula (TTJ) from Introduction ((2) from |17|). 
where II + = {II; x, X} symbolizes arbitrary "direct" 
process II together with conditions of its observation 
(external forces x(t) and thermodynamical forces X), 
II_ = {II; x, X} is time-reversed process in respect to 
11+ , AS = A5(II+) is system's entropy increment in 
the direct process, and X = eX . Of course, according 

to CEU), As*(n_) = -AS(n+) . 

C. From old to new relations 

It was already pointed out that left and right sides 
of equalities ([21]) and ([24]) describe evolutions of a sys- 
tem under consideration from one and the same equi- 
librium state, - more precisely, D eq (T,xo) (on the left) 
and D eq (T,xo) (on the right), - under mutually time- 
reversed protocols x(t) and ex(8 — t) . Obviously, if at 
that x(9) a;(0) then at least on one side the forces 
have a jump in the beginning of observation. 

This consequence of the reversibility is just what 
should be used if the system has no equilibrium states at 
x 7^ .To or has but we are interested in processes of transi- 
tion to such states from x — xq (see p.II.E). Otherwise, it 
is possible to replace (|2~Tj) and (|24|) by formally equivalent 
relations like "Crooks equalities" |34l - l38| where the men- 
tioned jumps disappear since these relations do compare 
evolutions from different equilibrium states, D eq (T, x(0)) 
and D eq (T,x(9)),. 

To write out such the relations, let us make in (f2"Tj) 
replacement 

${r(t)} => ${r(t)} x (25) 

x exp{-[H(T(9),x(0)) - H(T(0),x Q )}/T} , 

where <!>{r(i)} again is arbitrary functional. Then on 
the left in place of exp (-E/T) another exponential will 
appears, exp (— W/T) , while on the right the averaging 
over initial microstates will be performed with another 
weight function, 

exp {[F - H(T, x(0))]/T} = exp (-AF/T) D eq (T, ex{6)) 

(we took into account that Fq = F(xq) and F(ex) — 
F(x) ). Choosind here x(0) = xq , we obtain relation 

(${r(*)}exp(-W/T)) 1B(t) = (26) 
= exp(-AF/T) (a>{f(t)}) £(t) 



Now left and right sides correspond to protocols x(t) and 
x(t) (or, if indicating also initial states, {x(0),x(t)} and 
{x(0),x(t)}) which continuously join one with another. 
This is achieved thanks to difference between (equilib- 
rium) initial ensembles on the left and on the right. 
Next, insertion of (I22|) into (|26|) yields equality 

P[Q; x] exp {-W/T) = exp (-AF/T) P[Q; x] , (27) 

replacing (j2~4")l after the transition to new ensembles. This 
equality can be deduced also directly from (l24t with the 
help of multiplying (|24[) by the exponential factor from 
(f25*|) and exploiting the expansion (f22|) . 

D. More about interrelations of "new" and "old" 
results 

It should be emphasized once agin that the demon- 
strated transformation of of "old" relations (FDR) (sug- 
gested in 1977-1979) to "new" ones (suggested in 1997- 
1999) presumes that systems under consideration are 
closed, i.e. possess thermodynamically equilibrium states 
D eq (T,x) at any (constant) force values x. For open 
systems, - as we saw in p.II.E, - distributions D eq (T,x) 
do not exist (are non-normalizable) at x ^ . Conse- 
quently, all "new" relations, ([26]). (|27|) . etc., where W 
and/or AF appear, lose significance, like JE ^ does. 
At the same time, "old" relations do work with success. 

Especially it is necessary to dwell on the Marcovian re- 
alization of FDR. As it was shown in [TTJ] , generaled FDR 
are formally compatible with assumption about Marco- 
vian property of evolution and fluctuations of one or an- 
other set of variables Q(T(t)) , though really the "Mar- 
covianity" never takes place (if only Q does not coin- 
cide with r). This fact allows to formulate definite rules 
of construction of such Marcovian "stochastic models" 
what satisfy all FDR and therefore automatically agree 
with both reversibility of microdynamics and principles 
of (statistical) thermodynamics of irreversible processes 

0, d BEE il. 

In respect to closed systems, such rules in the main 
were considered alread y by Stratonovich (see [To| and ref- 
erences therein and in [ll], [l^, Gil ) who had started from 
the principle of detailed balance. We in [ll| extended his 
results to systems with time-dependent external forces. 

The corresponding theory (Stratonovich's and our) 
presented in 1977, if not earlier, completely include the 
Crooks theory which appeared first in Marcovian repre- 
sentation [34], |35| • Differences between these theories are 
in non-principal details only: terminology, mathematical 
design and notations. At that, transition from Marcovian 
language to general one formally adds nothing new, since 
deterministic (Hamiltonian) dynamics can be treated as 
limit case of Marcovian random process (with Q = T). 

In 1978-1981 we formulated generalization of the 
Stratonovich's Marcovian theory to open systems [H, 0, 
[i~9l ]. The corresponding our results, as far as we 
know, still are not reproduced by new authors. 
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Below, we shall speak about open systems. But before, 
since referring to Marcovian models, we have to discuss 
them at least briefly. 

E. Marcovian FDR 

Let a set of variables tp(t) = ip(T(t)) is sufficient to 
express through it the interaction energy (fTU| (for closed 
systems) and/or the dissipated power (for open systems), 
that is the integrand in (fT12j) . If ip(t) is assumed to be 
a Marcovian random process, then complete information 
about it is concentrated in kinetic equation for density 
P(t, tp) of its probability distribution: 

P(t,tP) = K(tP,V,x(t))P(t,tP) , (28) 

where V = d/dip , and K is differential or integral "ki- 
netic operator" (as, for instance, in the Fokker-Planck 
equation and Kolmogorov equation, respectively). For 
a closed system, this operator should satisfy operator- 
valued equality 

K(i/>, V, x] P eq (V»; x) = P eq x) K* {stp, eV, ex) (29) 

which accumulates all the consequences of FDR. Here, 
f is symbol of operator conjugation, or transposing, 
in the Sturm-Liouville sense, and P eq {ip',x) — J 5(ip — 
ip(T)) D eq (T, x) dT = P eq (etp;ex) is equilibrium distri- 
bution under constant forces. It represents stationary 
solution of (|28|) . 

If a system is open in respect to forces x , then (|2"9")) 
should be replaced by another operator-valued equality, 

K(^.y, x )p eq (ij) = 

= P e9 (^)[K t (eV,eV,ex),+N(VAx)/T] , (30) 

where P eq (%p) — P eq (tp,0) = P eq {£tp) is equilibrium 
distribution (existing at x = only), and N(t/>,x) is 
the dissipated power (x ■ I(tp) in the case (TH| - p5[) ). 
Now, stationary solution of ([28| at x ^ gives non- 
equilibrium distribution of tp . 

If there are simultaneously forces of both "closed type" 
and "open type" , then instead of equalities ((29)) and (|30|) 
one should write their obvious hybrid. 

Notice that is merely time-local form of the de- 
tailed balance (of direct and reversed processes of tp(t) 's 
interaction with thermostat), while (|30|) expresses vio- 
lation of the detailed balance to the extent of instant 
"entropy production" N(tp, x)/T . 

It is quite evident that equality (|30[) represents differ- 
ential form of the equality ([24]) (with Q = tp). Less 
evident fact is that equality (|2"9")l represents differential 
form of (|27|) . i.e. "Crooks equality". Nevertheless, it is 
so, and one can say that Crooks equality is expression 
of time- local detailed balance which (as shown in [TTJ) 
takes place even in strongly disturbed closed systems. 

It is necessary also to pay attention to degenerated 
Marcovian models where JC(ip, V, x) = — V • K(tp,x) . 



This means, obviously, that evolution of tp(t) is deter- 
ministic although may be dissipative, "with friction" (de- 
pendent on system's state), that is thermostat does not 
produce a noise ("has zero temperature"). At that, the 
distributions ("invariant measures") P eq (tp,x) become 
singular, while limit form of the relations ([29]) and (|30|) 
(in the zero noise limit) prescribes inversion of sign of the 
friction under time reversal. Therefore all FDR remain 
valid (at least in definite class of trajectories ip(t) ). Such 
models help to describe some numeric experiments [4l| 
which, however, are very artificial and far from actual 
problems of statistical mechanics. 

F. Fluctuation theorems 

Let us return to p.III.B and choose ${r(£)} = S(E — 
E{T(t),x(t)}) , where E{T(t), x(t)} is dissipated energy, 
(TT2")) or (TT5]) , considered as functional of phase trajectory 
of a system and applied forces (eventually, function of 
initial point T and functional of x(t)). Then average 
(f2"0"| turns to density P(E;x) of probability distribution 
of E , and relation (f2"Tj) gives us above mentioned equality 

09- 

The same follows from (|24)) and clear expression 
P(E;x) = J5(E - E[Q,x])P[Q;x}dQ, if notice that 
E[Q, x\ = — E[Q, x] . In quite similar way FDR (TTJ) im- 
plies analogue of (fT6| for distribution of AS : 

P(AS; x, X) exp (-AS*) = P(-AS; x, X) (31) 

Such particular FDR, transparently visible even with- 
out literal manipulations, in current literature are termed 
"fluctuation theorems" (FT) 0, . We in our time 

had written out and exploited them in somehow different 
forms, more practical at that time (see below). 

Further, we shall illustrate, by the examples from 
p.II.E, how FT were applied by us in 1979 [l2| and later 
to open systems. 

G. FDR for transport processes 

Consider, for demonstration of one of possible appli- 
cations of FDR, charge transport through a conductor 
under constant (after switching on at t = 0) voltage 
drop x . Hamiltonian of the system has form (TH| where 
Q(r) represents transported charge. We are interested 
in amount of transport during the observation time in- 
terval: AQ(9) = Q(T(9)) - Q(T) . The corresponding 
dissipated energy equals to E = xAQ . 

Let us combine exact FDR and a simple stochastic 
model of the system. The FDR will be delegated by 
relation (FT) (TTS]) written in the form 

P(AQ; x) exp (-xAQ/T) = P(-AQ; x) , (32) 

in fact used in jl2j |. with probability density P now re- 
lating to charge. What is for the model, assume that 
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our conductor is a contact (like e.g. p — n -junction) 
and therefore charge is transported through it by dis- 
crete portions ±e forming two opposite Poissonian ran- 
dom flows. This means that average value of the current , 
I(x) = (AQ)/0 , and spectral power density of the cur- 
rent noise, S(x) = [(AQ 2 ) — (AQ) 2 ]/9 , are expressed by 
formulae 

I = e[ra_|_ — n_] , S = e 2 [n + + nJ\ , (33) 

in which n± = n±(x) are mean numbers of the elemen- 
tary charge portions transferred per unit time in opposite 
directions. 

Clearly, FDR establish definite connection between 
n+(x) and n-(x) . In [12j | it was extracted from relations 
for characteristic function of AQ equivalent to (f3"2j). In- 
stead of this, here we merely can surmise that relation 
(f32|) is valid not only in respect to AQ as the whole but 
also in respect to elementary transfer events: 

n+(x) exp (— ex/T) = n-(x) (34) 

From this relation and (|33l) one obtains following relation 
between power of on-equilibrium noise and mean current 
(current- voltage characteristics) : 

S(x) = el(x) coth(ex/2T) (35) 

At e\x\ <C 2T it reduces to the Nyquist formula for "ther- 
mal noise" while in the opposite case to the formula for 
"shot noise". 

In such way FDR help to reveal universal connections 
between dissipative non-linearity of a transport process 
(I(x)), its noise characteristics (S(x)) and type of its 
statistics (here, simple Poissonian one). More compli- 
cated examples of this kind and corresponding more gen- 
eral FDR can be found in [IMl E3, . 

H. FDR and 1/f-noise 

Just considered stochastic model has principal defect: 
in it the elementary random events have a priori ("in 
advance") prescribed relative frequency (time-averaged 
number of events per unit time, or their "probability 
per unit time"), n±(x) , independent on concrete exper- 
iment, i.e. concrete phase trajectory of a system. Al- 
though, as it was shown by Krylov many years ago (42j, 
statistical mechanics gives no grounds for such assump- 
tions. 

One can understand this statement already on intu- 
itive level. Indeed, the mentioned assumption would be 
likely if the system remembered a number of past events 
and compensated its deviations from a "norm" by means 
of opposite deviation of number of next events. But this 
is impossible if the system forgets about events soon af- 
ter they happened. Then it does not distinguish between 
"norm" and "deviation" and therefore produces fluctua- 
tions in the number of events proportionally to its "nor- 
mal" (average) value. This means that relative frequency 



of events ("probability per unit time") undergoes low- 
frequency fluctuations with 1/f-type spectrum. 

For the first time similar reaso ning s were suggested 
and mathematically formulated in [43| and [2l|, [22[ and 
later confirmed on the base of statistical mechanics in 
[IIHIilllMIIii and other works, first of all in applica- 
tion to random walks ( "Brownian motion" ) of atomic-size 
particles. 

It should be underlined that the fluctuations ("1/f- 
noise") of relative frequencies by their very nature do 
not violate existing (anyway prevalent) balance or def- 
inite disbalance of mutually time-reversed events (pro- 
cesses). Therefore, - as it follows from generalized FDR 
[UHllIll, - various particular FDR like (O hold also 
for fluctuating relative frequencies and all derived "ki- 
netic" quantities. For example, the Einstein relation 
D = Tii between diffusivity and mobility of a walking 
particle, D and /j, , can be extended to their fluctuations 
[23| (as well as the Nyquist formula to fluctuations of 
conductance and fluctuations of "instant" spectral power 
density of "white" electric noise [U, HH). At that, to 
substantiate such statements, the relation (f3"2j) (formula 
(A4) from [23j) is quite sufficient. 

Due to these circumstances, - as it was suggested al- 
ready in [43j], - one can separate fast fluctuations (white 
noise) and low-frequency fluctuations (1/f-noise), making 
use of primitive phenomenological language but taking in 
mind its rigorous statistical- mechanical equivalent. Next, 
consider in such way statistics of random walk of probe 
("marked") gas particle, returning to the second example 
from p.II.E and basing on results of [H, HI, 0, |45| . 



I. Molecular Brownian motion and FDR 

Now, in the relation K32J) ((A4) from HH) AQ will de- 
note displacement, or path, of "Brownian particle" (BP). 
Let R be projection of AQ onto the external force di- 
rection. In the widely known simplest stochastic model 
of Brownian motion the particular FDR (FT) (l3"2j) is sat- 
isfied by the Gaussian distribution 

m^.P^^ fstSS^im ( 36) 

Here, it is assumed, of course, that the observation time 
9 is much greater than velocity relaxation time (mean 
free path time), r , of BP (probe gas atom). 

However, honest consideration of the exact BBGKY 
equations for infinite chain of many-particle distribution 
functions of a fluid shows that expression (|36[) is incom- 
patible with absence of (instant) statistical correlations 
between the BP and gas atoms far distanced from it. A 
true expression (first obtained in [44[ and then by differ- 
ent method in [28|, Hf|) can be represented by superposi- 
tion of the Gaussian distributions with various values of 
BP's mobility: 
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p CO 

P{R;x) = / P„(R;x)U e ( f i)dp 
Jo 



(37) 
(38) 



where S(-) is a "cut off" function which fast decreases 
at infinity and equals to unit at zero, S(0) = 1 , and Vo 
is characteristic thermal velocity of gas atoms (speed of 
sound). From here for variances of the BP's path and 
the dissipated energy E = xR it follows that 



(39) 



(R, R) = 2Tp 6 + {p x9) 2 F [\n - 
(E, E) = 2T(E) + (E) 2 F (in 



(40) 



Here and below the angle brackets with commas inside 
them (Malakhov's cumulant brackets) denote joint cumu- 
lant of random quantities separated by the commas. The 
second term in (|39p and (|40|) says about 1/f -fluctuations 
of mobility and dissipated power, respectively. At x = , 
similar asymptotic characterizes fourth-order cumulant 
of R , thus reflecting identical fluctuations of BP's diffu- 
sivity D = Tp [Illllll. 

The function Ug(p) is effective BP's mobility 

probability distribution. Its power-law long tail is gen- 
erally typical for distributions accompanying 1/f- noise 



At (E) 



fl Q X 



2 9 > T this cubic tail mani- 



fests itself in the path distribution (|3"?| on the right (if 
x > 0) : P(R;x) w (R) 2 /R 3 at R > (R) . Correspond- 
ingly, the same tail appears in distribution of dissipated 
energy in flU]) : P(E;x) w (E) 2 /E 3 at E > (E) . Hence, 
probabilities of "large deviations" of the path and dissi- 
pated energy highly exceed that predicted by the Gaus- 
sian model (|3l))) . 

Such distributions were repeatedly observed in experi- 
ments with non-stationary photo-currents (charge injec- 
tion currents) [46j |. 

It is interesting that most short way to these results 
runs from the FDR [28]. Let us choose in (HD 

${r(;)} = S(Q(9) - R) S(Q)H[1 + <f>(q k )} , 



where 4>{q) is some function of space coordinates, and 
the product is taken over all gas atoms (except the BP 
itself). Then, left side of (|21l) characterizes influence of 
initial spatial non-uniformity of gas onto BP's walk, while 
right-hand side describes statistical correlations between 
BP's path (during all the observation time) and current 
microstate of gas (in configurational space). Further, 
choosing the function <p{q) properly, one can extract 
from (|2"TT) relation 



= P(R;x) [is(p\R;x)-v]d 3 p, (41) 



dP{R; x) 
dv 



where v is mean density of gas, v(p\R\x) is conditional 
average value of gas density at distance p from BP under 
given its path, and for simplicity gas is assumed rarefied 
enough (besides, our designations here somehow differ 
from that in @]). 

From fiU) it follows that 



d In P{R;x)/d\n v > -i/fi 



(42) 



where f2 represents characteristic space volume to which 
the correlations of gas with BP do extend. On the other 
hand, in the Gaussian model (|36[) . in view of the known 
dependence p oc D cx 1 / v , we have 



d In P^{R;x) 
(9 In v 



AT 



(R)J 



Comparing this expression with inequality (|42p. we see 
that they certainly are incompatible if the "correlation 
volume" £1 is bounded above by a finite value. 

Hence, if the gas stays indiiferent to (forgets about) 
outcome, R, of BP's walk, then it is unable to supress 
large values of R so categorically as the law ([36]) does 
require. 

At the same time, the law (|37|) - ([38l) is well compatible 
with (14"2"1) . at fi = 2/v . Of course, this value of ft can 
not be obtained from FDR themselves only, its calcula- 
tion needs in the whole BBGKY hierarchy [23l [jij. 



J. Variance of dissipation fluctuations 

The previous paragraph gave example of large fluctua- 
tions of dissipation whose magnitude, according to (|4"0"j) . 
is of order of mean dissipation value. Another such exam- 
ple was considered in [24| . If energy is dissipated through 
not one but many, N 3> 1 , degrees of freedom, then vari- 
ance of dissipation fluctuations, along with magnitude of 
power- law tail of their distribution, will be approximately 
N times smaller. Anyway, for systems with Hamiltoni- 
ans of type (fT4"|) there are exact equalities [24| 

(E) = ^ [ x(l)(I(l),I(2)) x(tMt _ 2) x{2) dld2 , (43) 

1 Jl>2 

(E,E) = 2T(E) +| / x(l)x(2)x 

1 Jl>2>3 

x (7(1), 7(2), J(3)) B(t)l , ( t_3) z(3) dld2d3 (44) 

Here, the ciphers replace time arguments (and their num- 
bers), and T](t) is the Heavyside step function, that is in 
second- and third-order cumulants in the integrands the 
most early (right-hand) currents ( 7(2) and 7(3) , respec- 
tively) belong to still undisturbed (equilibrium) system. 
These equalities follow fro m g eneral FDR for cumulants 
[l2l ll4T ] (importantly, in [14| we had written out two 
variants of such FDR, for systems of "open" type and 
"closed" type separately). 
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In many applications one can suppose that the third- 
order cumulant (or result of its integration) in (l4"4l van- 
ishes at x = . Then, under weak disturbation the first 
term in (|44p is oc x 2 while the second oc , and 

(E,E)=2T(E) +1 / x(l)x(2) x 

1 J2>3,4 

x (7(1), G(2, 3), 7(4))o x(3) x{A) dld2dMA , (45) 

where G(2,3) = [6I(2)/Sx(3)] x =o is dynamical (ran- 
dom from the ensemble viewpoint) linear response of the 
currents to forces (24|. However, the second term does 
not reduce to a small correction of the first, if it grows 
with observation time faster then linearly. For instance, 
when it grows, like in (l4"0"|) . approximately oc 2 . At 
that, the integral in (|45|) determines values and spectra 
of equilibrium 1/f -fluctuations in kinetic coefficients of 
the system and power dissipated by it. 

IV. CONCLUSION 

We have presented to readers a brief review of the "gen- 
eralized fluctuation-dissipation relations (theorems)", for 
the first time obtained in 1977 in our work (Tlj . in com- 
parison with much later analogous results, in particu- 
lar, the "fluctuation theorems". Both our "old" and 
later "new" results express important consequences of 
time symmetry (reversibility) of microscopic motion, sub- 
stantially determining statistical and dissipative proper- 
ties of thermodynamically non-equilibrium physical sys- 
tems. However, the expounded comparison gives suffi- 



cient grounds to state, - in contrast to the misunderstand- 
ings observed in current literature (see Introduction), - 
that "new" results do not contain a principal novelty and 
in the sense of generality are significantly inferior to "old" 
results (if not per se reproducing them or even earlier re- 
sults by Stratonovich) . 

The matter is that our approach (reflected in the 
present paper), - based on constructive play with sta- 
tistical ensembles and Hamiltonians, - is more flexible 
and fruiful. Due to it, our "generalized fluctuation- 
dissipation relations" (FDR) from the very beginning 
were represented in the form foreseeing their applications 
to systems of not closed type only but also open type, 
while "new" results were oriented to the first of these 
two and lose meaning in respect to the second. 

In the framework of our "old approach" any "new re- 
sult" can be obtained from the "old" FDR, but the op- 
posite in the framework of "new approach" is in general 
impossible. Nevertheless, some of new formulations of 
FDR indeed are more comfortable in their domains of 
applicability and thus enlarge collection of practical ap- 
plications of FDR. 

On the whole, the generalized FDR bring all neces- 
sary tools for construction of thermodynamically correct 
stochastic models of real processes and systems. Regard- 
less of degree of complexity or roughness of a model, ob- 
servance of FDR at its level ensures its qualitative agree- 
ment with rigorous statements of statistical mechanics 
(and sometimes even leads closely to quantitative agree- 
ment, as was demonstrated, in particular, by examples in 
last paragraphs of the present paper) . Hardly this useful 
potential of FGR some day will be exhausted. 
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